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Abstract 

We apply the nonlocal condensate formalism to construct generalized sum rules (including 
0(a s )-radiative corrections) for 7T-, p- and p'-meson wave functions of twist 2. Besides, we predict 
the lepton decay constant f p > and estimate the mass of //-meson. For all these mesons we obtain 
the first 10 moments of longitudinal wave functions, suggest the models for them and discuss their 
properties. We consider the peculiarities of the QCD sum rules with nonlocal condensates for 
transverse p- and &i-meson wave functions. These results are compared with those found by Ball 
and Braun and by Chernyak and Zhitnitsky. 
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1 Introduction 



An important problem in the theory of strong interactions is to calculate hadronic wave functions 
<Pn(x), (fp(x), <pn(xx, X2, £3), etc. from the first principles of QCD. These phenomenological distri- 
butions of partons on the fraction xP of a hadron momentum P are a natural result of "factorization 
theorems" applied to hard exclusive processes [0, |, g ||. They accumulate all the necessary infor- 
mation about non-perturbative long-distance dynamics of partons in hadrons. 

In the standard QCD sum rule (SR) calculation of light meson wave functions (WF's), first 
introduced by Chernyak and Zhitnitsky (C&Z) [|| and recently re-estimated by Ball and Braun 
(B&B) for the p- meson ||, it is assumed that the correlation length A of vacuum fluctuations is large 
compared to a typical hadronic scale ~ l/m p . Thus, one can replace the original nonlocal objects 
like M(z 2 ) = (q(0)E(Q, z)q(z)) [] by constant quantities of (?(0)g(0))-type. Based on this hypothesis, 
the well-known QCD SR approach j7| has been applied in || to calculate the first two moments 
= Jq (p(x)(2x — l) N dx with N = 2, 4 for WF's of light mesons. And just from these moments 
the whole WF's have been reconstructed which are now referred to as C&Z WF's. 

Now it is known that hadronic WF's are rather sensitive to the widths of the function M(z 2 ) 
and of other nonlocal condensates j|, ^, [h| and the crucial parameter A • m p ~ 1. Therefore, one 
should use nonlocal condensates (NLC's) like M(z 2 ) whose forms reflect the complicated structure 
of QCD vacuum. Certainly, these objects can subsequently be expanded over the local condensates 
(q(0)q(0)), (<?(0)V 2 g(0)), etc. (here V M = d p — igA p is the covariant derivative), and one can come 
back to the standard SR by truncating this series. Our strategy is to avoid an expansion of that sort 
because we thus lose an important physical property of non-perturbative vacuum - the possibility of 
vacuum quarks and gluons to flow through vacuum with non-zero virtuality k 2 ^ 0. Indeed, the 
average virtuality of vacuum quarks (k 2 ) = X 2 is not small whereas the standard approach inevitably 
suggests that k is not small but exact zero. The value of X 2 can be extracted from the QCD SR 
analysis 0] and is connected with the condensate of the next (d = 5) dimension 

m 2(55) k 

that is of an order of the typical hadronic scale, m 2 ~ 0.6 GeV 2 . An estimate for X 2 in the framework 



of the instanton liquid model yields a similar number, and the lattice QCD calculations [13] give 
A 2 = 0.55 ± 0.05 GeV 2 , a similar estimate A 2 = 0.5 ± 0.05 GeV 2 has been obtained injpj. 

Careful inspection of consequences of that approach to QCD SR for pion WF []|, ||, |15| has revealed 
that the introduction of the correlation length A ~ 1/ X q into condensate distributions produces much 
smaller values for the first moments of pion WF than C&Z values. As a result, the pion WF is 
strongly different in shape from the C&Z and approaches the asymptotic WF (p as (x) = 6x(l — x) pj, 
i.e. , ^(x) « ip as (x) Later, this WF was confirmed by independent consideration of the QCD SR 
directly for the function ip n (x) based on the non-diagonal correlator and on the advanced smooth 



distribution function for the quark nonlocal condensate [18, 19]. The key element of both the ways 
was to take into account the main physical reason - vacuum correlation length A ~ 1/A 9 is of an 
order of l/m p . 

Our goal here is to show that in the case of QCD SR for the p-meson channel the situation is similar 
to the pion one. All the predictions of the standard QCD SR (the C&Z ones for the longitudinal case 
H and B&B ones for the transverse case || - we call them the Local QCD SR) for moments (£, N ) P 
with N > 2 could not be considered as reliable; as we show in our analysis, even the value of p is 
indefinite. We apply the NLC formalism to calculate the diagonal correlators for />meson currents, 
introduced in ||, and construct generalized SR including 0(a s )-radiative corrections to obtain WF's 
of twist 2. The first ten moments of longitudinal WF's of p- and //-mesons are estimated (Table 1) 
and the models for them (see FigsJ^-Q) are suggested. Also, we predict the lepton decay constant 
fpi and estimate the mass of //-meson. We discuss the strong dependence of results for moments of 

1 Here E(0,z) — Pexp(i dt^A^ifyTa) is the Schwinger phase factor required for gauge invariance. 
2 Note here other independent evidences that the pion WF is close to its asymptotic form: lattice calculations of the 
corresponding (£ 2 },r in |l(|; QCD SR estimate of the magnitude of at the midpoint x = 1/2 in Jl7| . 



transverse WF's in p-channel on the model of nonlocal gluon condensate contribution. The calculation 
technique is the same as in Refs. ||, |9|; therefore, the corresponding details are omitted below, but 
we shall through the text keep the connection with the well-understood pion case. 



2 Generalized sum rules for the tt- and p-meson channels vs the 
standard version 

For the helicity-zero charged vector meson M^\=o = Mi, Ml = p, p' the leading-twist WF is 
defined as 

(0 | d(zh,u(0) | M L (p))\ z2=Q = ffarfidxe*™ ^/ L (x) (2) 

(In accordance with the B&B definition, in this case e^ =0 (p) — Viil m M h as p z —>■ oo, where e fl (p) is 
the polarization vector); and for axial mesons (M A = tt,A\ . . .) as 

(0 | d(z) lfll5 u(0) | M A (jp)) = if MAPtl [ dxe ix W <pm a {x). (3) 

2 =0 Jo 

The information about <p L {x) can be obtained from the correlator I( n o)(q 2 ) of vector currents V( n )(y), 
see e.g. [§: 

i Jdye^(0\T{v+ ) (y)V (n) (0)}\0)=z n (z q T +2 I {n0h (4) 

and the information about ip (x), from an analogous formula with the substitution V — ► A and 
Ar n \(y) = d{y)z^ (zV) n u(y). The corresponding formula for the unit helicity state, M± = Mm =1 , 
M± = p, b\ . . . , is as follows 

(0 | d(z)a^u(0) I M ±(P))\ z2=0 = */mx ( e »^ ~ £ uip)P„) f Q dxe ix ^ tp T M± {x). (5) 
Moments of these functions tp T (x) are extracted from the correlator J( n o)(q 2 ) of tensor currents 

i J dye^(0\T {T$(y)T£ n) (0)} |0> = -2i n (zq) n+2 J (n0) , Tfc(y) = d(y)a^z a (zV) n u(y). (6) 

This correlator J( n o)(q 2 ) contains the contribution from states with different parity (see the analysis 
in ||). Therefore, the contamination from 6i-meson (^J PC = l" 1 ^ in the phenomenological part of 
SR is mandatory. 

The theoretical "condensate" parts for both the correlators, J( n0 ) in (Q) and J( n o) in (g), contain 
the same 6 universal elements as for the pion case. An exception concerning transverse sum rules 
(see (B)) will be discussed below. The diagram origins of these elements A$r (x',M 2 ), where M 2 
is the Borel parameter, are described in detail in ||, [|. By the direct SR formulation for WF, one 
immediately constructs a "daughter SR" for any functional of <pm( x ) ( n °t only for moments (£ )). 
Let us write down the final SR's including WF's of p-meson and next resonances p' and b\ into 
the phenomenological parts (the corresponding SR for the pion (axial) channel is written below for 
comparison): 

{U) 2 ^{x) + {f Al ) 2 ^ Al (x)e- m -J M2 = f° P T\x-s)e-°l M2 ds+ (7) 

Jo 

+ A$ G (x;M 2 ) + A%(x;M 2 ) + A$ y (x;M 2 ) + A$ Tl (x;M 2 ) + A$ r2 (x;M 2 ) + A$t 3 (x;M 2 ) ; 

[ff)\ L M^ IM2 + {fi)\ L A*)e- m *' /M2 = [J PT(*; s)e-^ 2 ds + (8) 



+ A$ G (x;M 2 ) - A$ s (x;M 2 ) + A<Z> v (x;M 2 ) + A$ Tl (x;M 2 ) + A$ Ta (x;M 2 ) + A<S> n (x; M 2 



(/J) V(*)e _m ' /Ma + (/£) ^K"^ = / ° s)e-l*ds + (9) 

+ A$ G (x;M 2 ) - A$' G {x; M 2 ) + A$ v (x;M 2 ) + A$ Tl (x;M 2 ) + A$t 2 (z;M 2 ) - A<5>t 3 (x; M 2 ) , 



where Sq' L ' T are the effective continuum thresholds in the pion and the p-meson L and T cases, 
respectively. Perturbative spectral densities (f^ r T (x; s) are presented in an order of 0(a s ) in g § 
for the L case and in || for the T case (see Appendix B). Radiative corrections reach 10 % of 
ff^ris ~ 1 GeV 2 ). Contributions A$ r (x;M 2 ) depend on a specific form of NLC's M(z 2 ), etc. 
To construct SR for WF's, it is useful to parametrize these NLC behaviors by the "distribution 
functions" [||, |9], [L8], a'la a-representation of propagators, e.g. , fs(v) for the scalar condensate 
M(z 2 )0 

P OO P OO /^OO 

M (z 2 ) = (g(O)g(O)) jf e™ 2 /* f s (u) dv, where jf f s (u)du = l,j uf s (u)du = -f. (10) 

The function fs(v) and other similar functions /r(^) describe distributions of vacuum fields in virtu- 
ally for every type of NLC. They completely determine the r.h.s. of SR's in (|7])-(§). The general 
forms of elements A<I>r(x; M 2 ) as functionals of fr(v) will be published in a separate paper. For the 
standard (constant) condensates (G(O)G(O)) and (q(0)q(0)} these distributions are trivial in form, 
e.g. , fs{v) = 5(y) (Appendix A). To include the condensates of the next (~ A 2 ) and higher dimen- 
sions into consideration , one should add the contributions to fs{ v ) proportional to the derivatives of 
5-functions, 5{y) , 5(y) , . . .. It is clear that since there is no QCD vacuum theory, merely models of 
real distributions can be suggested for these distribution functions. However, for the purpose of QCD 
SR's for moments (£ N ) we need an approximate information about the fviy) behaviour. Therefore, 
we here apply the simplest ansatz || ||, like fsiy) = S(v — A 2 /2), to take into account only the main 
effect, the non-zero average virtuality (A; 2 ) of vacuum fields. This form of fsiy) can be regarded as a 
result of resummation of an infinite subset of the above-mentioned contributions ~ (A 2 /2) n 5{y)^ 
connected with the single scale A 2 [||, [|. The corresponding expressions for A^G,s,v,Ti{x; M 2 ) are 
collected in Appendix A. 

Now let us take the limits A 2 -> 0, A<S> T {x,M 2 ) -> Aip v (x,M 2 ) and p^ rt (x,s) -> pf orn (x,s) 
for SR in eq.(||) to return to the standard approach (see these reduced elements in Appendix A). 
We try to inspect the subtle points and the range of validity of C&Z SR. These authors extracted 
(£ ) exactly in the same way as the f p value (B&B limited themselves only to extraction of (£ 2 )). 
However, the nonperturbative terms in their sum rule (the /(/-contribution is omitted for simplicity) 
have a completely different A-dependence compared to the perturbative one and, a priori, it is not 
clear whether a straightforward use of the 'A = technology' can be justified for higher A (for 
definiteness, we consider here only the /9-meson (longitudinal) case; the same arguments apply also 
to the pion case, see criticism in || 10, [D|) 



\ Jp ) ^ lp ^ 4vr 2 (A + l)(A + 3) 

3M2 + %^ + ^-7)<«. (ID 



4vr 2 (A + l)(A + 3) 12M 2 81 v M 4 

The scale determining the magnitude of all hadronic parameters including sn (the "continuum thresh- 
old" Qj) is eventually settled by the ratios of condensate contributions to the perturbative term. If 
the condensate contributions in the C&Z sum rule ( |TTj ) had the same A-behavior as the perturbative 
term, the A-dependence of (£ N ) would be determined by the overall factor 3/(A + 1)(A + 3) and 
the resulting WF <p(x) would coincide with the "asymptotic" form ip as {x). 

However, the ratios of the (qq)- and (GG)-corrections to the perturbative term in eq. (|Tl|) are 
growing functions of A. This reduces the predictable power of the local QCD SR's with the growth of 
A. To reveal consequences of this effect more clearly, let us consider the so-called SR fidelity windows, 

3 In deriving these sum rules we can always make a Wick rotation, i.e., we assume that all coordinates are Euclidean, 



i.e., regions of the Borel parameter M 2 where one should obtain reliable SR predictions. In accordance 
with the QCD SR practice 0, these fidelity windows are determined by two conditions: the lower 
bound M 2 by demanding that the relative value of (GG)- and (qq) -contributions to OPE series 
should not be larger than 30%, the upper one M 2 . by requiring that a relative contribution of higher 
states in the phenomenological part of SR should not be larger than 30%. Suggesting independence 
of the threshold of N (s N ps s ss 1.5 GeV 2 ), we have in the case of N = 0: M 2 = 0.4 GeV 2 , 
M 2 = 1.34 GeV 2 . But for N = 2 we have M 2 = 0.73 GeV 2 , M 2 = 1.34 GeV 2 , and for N = 4 - even 
M 2 = 1.5 GeV 2 , M\ = 1.34 GeV 2 . That is, the fidelity window shrinks to an empty set in the last 
case. C&Z had to suggest that S2 ~ 1.9 GeV 2 and 54 ~ 2.2 GeV 2 to extend the fidelity windows and 
to obtain any stability with respect to M 2 . It is difficult to imagine such a strange type of a spectral 
model, but there are no principal objections. 

The situation with (£ 2 )p is to a certain extent special: though there exists a non-empty fidelity 
window of Local QCD SR for this quantity, the stability of SR in this window is rather poor and 
this is due to overestimation of vacuum condensate contributions. To improve the situation, C&Z 
suggest that the continuum threshold should be increased, which means assigning a part of higher 
state contributions to the ground state contribution, that is, overestimation of (£ 2 ) p . On the contrary, 
the authors of Q put the threshold unchanged, S2 = sq, and give the prediction without establishing 
any stability as a mean value of SR at M 2 = 0.8 GeV 2 and at M 2 = 2.0 GeV 2 ; the overestimation 
of (£ 2 } p in this case is due to that of the condensate part of SR. 

In our opinion, there is no need to propose such an exotic spectral model (sn = so + const • N) 
because the reason for this "exploding" behaviour of Local SR is quite evident, namely, a completely 
different dependence on N of the perturbative (the first term in the second line of Eq.(|TT|)) contri- 
bution and of condensate ones. And the origin of this difference was explained in a series of papers 
H 10, 15]; this is due to the Taylor expansion of initial nonlocal objects like (q(0)E(0, z)q(z)) in 



powers of z. The first constant term of this expansion, (qq), produces an (iV)°-dependent term in 
SR ([ll]); the next term, an (A^) 1 -dependent, and so on. 

On the contrary, the NLC terms A&r(x; M 2 ) in (||) and (||) lead to the moments (£ N ) which well 
decay with iV growing; so a physically motivated ^-independent continuum threshold Sq naturally 
appears in the SR processing. 



3 The moments and the models of wave functions 

Before analyzing the results of processing of SR (|8|) and (||) for the moments ($, n )m, let us consider 
the peculiarity of the QCD SR structure represented in the RHS's of Eqs. ([?]-||). Unlike the 7r-meson 
case, the contribution of the numerically most significant "four-quark condensate" A<]?s (x,M 2 ) Q 
is equal to zero (for the T case, see also the sign of A$>t3) or even has the opposite sign (for the 
L case). For this reason, the role of the vacuum interaction for the p-meson is weaker than for the 
pion. As a consequence of such an SR structure the nonlocal effects partially compensate themselves. 
Therefore, the extracted values of (£ 2 )a/ in the framework of NLC SR don't drastically differ from 
the results of B&B obtained in the standard way 0. However, the sensitivity and stability of NLC SR 
are much better than for the standard one, compare the accuracy for the first moments in Table 1. 
This allows us to estimate the first ten moments in the L case for p- and //-mesons. For clarity 
see Fig.|l|, where the curves for (£ 2 }p are represented as functions of the Borel parameter M 2 in the 
range of stability window for two types of QCD SR's: the solid line is determined from NLC SR 
(||); and the dashed one, from B&B SR || (actually B&B paper contains only the graphics for the 
coefficient , but it is straightforward to obtain the corresponding curve for (£ 2 )p = 0.2 + ; 
for comparison we also depict in Fig.[l] the short-dashed curve for the NLC SR prediction without 
p'- meson contribution, because B&B SR has been processed just in this manner). As one can see, in 
the local SR case there is no stability at all. 

We have determined the following values for practically TV-independent continuum thresholds: 
s(j ^ 2.4 GeV 2 . Fidelity windows for the L case are: 0.6 GeV 2 < M 2 < 2.0 GeV 2 for all N = 

4 Note here that the results of the T case in an original CZ work should not be taken as a pattern for a standard 
SR, there is an error in the sign of the quark condensate contribution, see pq, 



0, 2, . . . , 10. The ranges of stability within these fidelity windows almost coincide with these windows, 
starting for all A at M 2 ps 0.8 GeV 2 . 







The moments (£, N )m(P' 2 ) at y? ~ 
(errors are depicted in brackets following a 


- 1 GeV 2 

standard manner) 


Type of SR 


f M ( GeV 2 ") 


A = 2 


A = 4 


A = 6 


A = 8 


TV = 10 


Asympt 


WF 


1 


0.2 


0.086 


0.047 


0.030 


0.021 


NLC : 


7T 


0.131(2) 


0.25(1) 


0.110(7) 


0.054(3) 


0.031(2) 


0.0217(7) 


C&Z : 


7T 


0.131 


0.40 


0.24 


- 


- 


- 


NLC : 


P L 


0.201(5) 


0.227(7) 


0.095(5) 


0.051(4) 


0.030(2) 


0.020(5) 


B&B : 


P L 


0.205(10) 


0.26(4) 


- 


- 


- 


- 


C&Z : 


P L 


0.194 


0.26 


0.15 


- 


- 


- 


NLC : 


P' L 


0.175(10) 


0.226(10) 


0.145(7) 


0.106(5) 


0.082(4) 


0.064(4) 


NLC : 


P T 


0.169(5) 


0.325(10) 


? 


7 


? 


? 


B&B : 


T 

P 


0.160(10) 


0.27(4) 










C&Z : 


P T 


0.200 


0.15 


< 0.06 








NLC : 


bj 


0.181(5) 


0.090(5) 


? 


7 


? 


? 


B&B : 


bl 


0.180-0.170 













Table 1. The moments (£, N )m{p- 2 ) at // 2 ~ 1 GeV 2 . C&Z give all moments normalized to 
the normalization point hq = 500 MeV; here we present these moments normalized to the 
normalization point /j, = 1 GeV. 

For transverse SR @, in its theoretical side, we face the problem that an important part of 
nonlocal gluon contribution (x, M 2 ) is not yet estimated. We have suggested a naive model, 

2 n (a s GG) 
) 6vrM 2 

2A_ A J („ ^\0(A<x)e{x<l-A) 



instead of the local result Aip' G yx, M j = q we P u * 



A& G fx, M 2 ) = A<p' G fx, M' 



1 - 2A 



This simulation eliminates end-point (x = 0,1) effects due to the influence of the vacuum gluon 
nonlocality, which is inspired by the analysis in [^] and our experience in the nonlocal quark case 
(see Appendix A). For the zero-th moment SR (A = 0) any model is not too important, they give 
a gluonic contribution close to one in the local SR, and we obtain the best stability with Sq ~ 2.3 
GeV 2 and for 0.75 GeV 2 < M 2 < 1.9 GeV 2 . But for the next moment the results are very sensitive 
to the model of gluon contribution. Indeed, these results (see Table 1) drastically change if we leave 
the local expression Aip' G (x, M 2 ^j unchanged: {£ 2 ) T p = 0.231(8) and (£ 2 )£ = 0.220(8)). So, we 
realize that our modelling is of great importance for the results discussed and we can't give a reliable 
prediction for (£, 2N )'^f )1 (A = 1, . . .) without essential calculational efforts. 



We also re-estimate the first ten moments for the pion case and confirm the main previous 
results of nonlocal approach ||. In this axial channel (see (R)) we include the pion and ^i-meson 
into the phenomenological part of SR and obtain: Sq ps 2.3 GeV 2 , and the ranges of stability are 
0.75 GeV 2 < M 2 < 1.9 GeV 2 . 

Another evidence of the efficiency of NLC SR is the estimate of the p'-meson mass. First, the 
//-resonance with tabular mass m p > = 1465 MeV was inserted into SR (|8|) to improve the stability. 
But at the second step m p > was estimated from our SR, see Fig.§. It appears to be rather close to 



the experimental value [22 



m theor = ± 37 Me y ^ m exp = ± ^ MgV _ ^ 

Possible models of WF's corresponding to the moments in Table 1 have the form 

tfTW 2 ) = / s (x)6.739 (/ s (x)) L305 (x 2 + (l-x) 2 ) 4 - 01 , (13) 

rt mod {x^ 2 ) = <f as (x) fl + 0.077 -Cl' 2 (0- 0.077 • C 4 3/2 (0) , (14) 



L,mod, 



*,/x 2 ) = ^ s (x) (l + 0.075 -C 2 3/2 (0+ 0.4 -Cl /2 (0- 0-04 -C 6 3/2 (0), (15) 



where £ = 1 — 2x, C^(£) are the Gegenbauer polynomials and fi 2 ~ 1 GeV 2 corresponds to an average 
value of M 2 . The important problem of sensitivity of the results to the condensate parameters in 
SR's will be addressed to in a separate paper. Here we only note that the main qualitative conclusions 
concerning the shapes of fppi and <p™ do not change under accessible variation of A 2 . 

To check the reliability of these models, let us estimate the functional /fojvf] = / dx 

Jo x 

that is often used in calculations of different form factors. This integral naturally appears in the 
perturbative QCD approach to the electromagnetic form factor F n [jl], Q], to the evaluation of 77* — ► 
7T |23j and also to the QCD description of heavy meson semileptonic decays |^4j . It is clear that 
/[</?m] is a new independent (of moments (C N }m) quantity. Besides, the values of I[^Pm] allow us to 
better discriminate between different models for the same <pm( x )- The I[v?m] can be obtained in two 
different ways, 

• from QCD SR adapted to the quantity I[<p L ] by integration with weight 1/x; 

• by direct integration of the WF models (|l3|)-(15). 



Let us first consider the results for the pion case. The range of stability here is as wide as 0.6 GeV 2 < 
M 2 < 1.9 GeV 2 and the quality of stability is very high. As it is seen from Table 2 (columns 3 and 
4), our estimates for /[i^-], obtained in both the mentioned ways, are consistent with each other 
within 10% accuracy. Moreover, the estimation for I[ip w ] from SR, is quite close to the one obtained 
independently by using the nondiagonal NLC QCD SR in |l^] (column 5). So, we may conclude 
that the nonlocal approach leads to self-consistent results which are not far from the estimate of 
I [fir] ~ 2.4 in [23] [] (see column 6) and can be reliably discriminated from the C&Z ones. 



A similar situation holds also for the /^-results, the stability is also of high quality in the range 
0.5 GeV 2 < M 2 < 2.0 GeV 2 . Agreement of both kinds of the results for I[<p L ] is rather good for the 
Pi-case (the discrepancy is smaller than 5%) and worse for the p L -c&se. The models (fl^)— (|T^) confirm 
the property of NLC SR's that WF's of a meson in the "ground state" approaches the asymptotic 



WF (due to nonlocality effects), proposed by Radyushkin [15]. The curve ^ m od( x ) ( see FigJH) is 



not far from the asymptotic WF curve and looks close to the naive B&B model. Note, nevertheless, 
that the latter visual closeness seems rather crude and does not allow quantitative conclusions, e.g. , 
NLC SR provides the values of moments (£ 2 }p, / [ l f L \ in 15% smaller than in the B&B case. One may 

5 Another model for tp w , ip™ odl — (p at '(x)(l + -kjC\ 12 '({)), corresponding to the first three moments, has been suggested 



in I 

6 



The authors of [E3[ state that the value I ~ 2.4 may be underestimated not more than by 20% 



expect that WF's of resonances would oscillate by analogy with the pion resonance tt' case fl8| , 19 1. 



Indeed, the shape of WF ip^, ,mod (x) (see FigH) for resonance looks similar to WF of tt'. 





WF 


SR 


WF 


WF 


SR 


WF 


WF 




(asymp) 


(here) 


(here) 


(©) 


(©) 


(B&B) 


(C&Z) 




3 


2.75 ±0.05 


3.03 


2.8 


2.4(< 2.9) 


— 


5.0 




3 


3.1 ±0.1 


3.0 






3.54 


4.38 




3 


4.7 ±0.2 


4.3 











Table 2. On the analysis of SR self-consistency, based on evaluations of the func- 
tional I[<pm] = / — dx in different approaches. 



o x 



4 Conclusion 

Our basic interest in the present paper is to explore the well-working method of NLC QCD SR in 
analyzing WF's in the /3-meson vector and tensor channels. As has been noted in the previous papers 
|9], |n], [l^, 18 , [l{| , just in problems of "nonlocal" characteristics of hadrons, such as wave functions, 



form factors, etc. , one should use the formalism of nonlocal condensates. Let us summarize the main 
results of this paper: 

1. The generalized sum rules for WF's of the p-meson and related resonances with nonlocal con- 
densates are constructed. Using the simplest ansatz Q for nonlocal quark condensates we 
obtain new estimates for the first ten moments of the p-meson and its resonance WF's. It 
should be emphasized that analogous evaluation within the standard QCD SR approach is 
impossible. 

2. We suggest the models for (see Figs. 3-4) longitudinal WF's of tt-, p- and //-mesons and verify 
their self-consistency. The form of the obtained longitudinal p-meson WF is not far from the 
asymptotic WF (this conclusion noticeably differs from the results of the Local QCD SR ||). 

3. As a by-product, we predict the lepton decay constant f p / = 0.175 ± 0.005 GeV 2 and estimate 
the mass of the p'-meson, m„i = 1496 ± 37 MeV, which is now under experimental investigation 
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4. For 7T- and /^-channels we perform processing of SR directly for the integral I[(p] = [ ^ M dx 

Jo x 

and obtain high-stable predictions for these quantities. Their values unambiguously favor WF's 
close to the asymptotic form rather than to the C&Z-one. Moreover, the result for 7r-meson 
is also consistent with the model WF obtained in a special variant of NLC QCD SR's (in 



the non-diagonal QCD SR approach fllq , IS] which directly connects 7T-WF only with nonlocal 
condensates) and with analysis of the 77* — > tt transition form factor . 



So we can conclude that the nonlocal QCD SR approach is self-consistent and gives the results of 
high stability. To complete the calculation for transverse p-meson WF, we need an exact result 
for the nonlocal gluon contribution to SR. An open problem of this approach is the determination 



of well-established models of distribution functions /r(^) from the theory of nonperturbative QCD 
vacuum. 
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Appendix 

A Expressions for nonlocal contributions to SR 

For vacuum distribution functions /r(^) we use the set of the simplest ansatzes 

fsiy) = 5{y-\ 2 q /2) ; f v (u) = 5' (v - \ 2 q /2) ; 
/ Ti (ai,a 2! «3) = 5 (qi - \ 2 q /2) 5 (a 2 - \ 2 q /2) 5 (a 3 - X 2 g /2) . 



(A.l) 
(A.2) 



Their meaning and relation to initial NLC's have been discussed in detail in [Q, |j| The contributions of 
NLC's A^r{x, M 2 ) corresponding to these ansatzes are shown below; the limit of these expressions to 
the standard (local) contributions ipr{x,M 2 ) - A 2 — > 0, A^r(x, M 2 ) — ► A(pr(x, M 2 ) are also written 
for comparison. Here and in what follows A = A 2 /(2M 2 ), A = 1 — A: 



A$ 5 (x,M 2 ) 



As 18 
M 4 AA 2 



{0 (x > A > x) x [x + (A - x) In (x)\ + (x — > x) + 
+9(1 > A)0 (A > x > A) [A + (A - 2xx) ln(A)] } 



A(p s [x, M 
A$y (x,M 2 
Aip v (x,M 2 
A$ Tl (x, M 2 

A<p Tl (x,M 2 
A$ T2 (x, M 2 

A(p T2 (x, M 2 
A$t 3 (x, M 2 



^9 (*(*) + (*-*)); 

^(x5'(x-A) + (x^x)) 
As 

{x5 (x) + (x -» x)) ; 

_3As 
M 4 



|[<5(x - 2A) - 5(a; - A)] - 2^j 9(1 > 2A) + 9(2A > 



9(x > A)9(x >3A-F 



A 



3x 1 + x 

A" ~ 6 A~ 



+ (x — > x) , 



3A S 
M 4 



(5' (x) + (x^ x)) ; 



_ |5(x_2A) 0(i > _ ^ (2A > x ^ > > 3A _ ^ 



1 + 2x - 4A 
AA 2 

2A S 



+ (x —* x) , 



M 4 
3A s x 



M 4 AA 

+ (x — ► x) , 



(x5' (x) + (x — > x)) ; 

9(2A > x)6(x > A)9(x > 3A - 1) 



x A 
2 ~A~A 



(A.3) 

(A.4) 
(A.5) 

(A.6) 



(A.7) 



(A.* 



A^r 3 

Aip G 



G 



G 



x,M' 



x,M' 



x,M' 



x,M' 



x,M' 



8vr 



3A S 
M 4 

(a s GG) 
24vrM 2 
(a s GG) 
24vrM 2 

(a s GG) 6 (A < x)#(x < 1 - A) 



(5 (x) + (x — > x)) ; 

{5(x - A) + (x -> x)) , 
(5 (x) + (x — >■ x)) ; 



6vrM 2 
6vrM 2 ' 



1 - 2A 



(A. 



(A.10) 



Here As = — {y/a^q(0)q(0)) 2 , for quark and gluon condensates we use the standard estimates 
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<v^?(o)?(o)> « (-°- 238 GeV ) 3 > (asGG) 



12vr 



0.001 GeV 4 @ and 



2 _ (q(ig^uGn^ 



0.4 ± 



0.1 GeV 2 , normalized at jj? « 1 GeV 2 . 



B Expressions for perturbative spectral densities 

First, p(x,s)^ rt in an order of 0(a s ) was calculated in []|, but there was omitted the trivial 

-term in correlators here we restore it. The 



term 2 In 



that follows from the In 



corresponding term for the T case p(x, s)^ rt has recently been presented in ||. We have recalculated 
it and confirmed this result. 



2vr 2 



XX. 



p P L r \x,s) - -—j.r.r [ 1 + 

pert / \ 



2vr 2 " 

3 f 

— jXX < 1 + 

27T 2 47T 



a s {p?)C F 

47T 

a s (fi 2 )C F 



2 In 



2 In 



7T 



+ 5 h ln 2 (x/x) 

3 



7T 



+ 6 — + ln 2 (x/x) + ln(xx) 



(B.l) 
(B.2) 

(B.3) 



Here /u 2 ~ 1 GeV 2 corresponds to the average value of the Borel parameter M 2 in the stability 
window; a s ( lGeV 2 ) ra 0.52. 
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Fig. 1: Curves (£ 2 )p as functions of Borel parameter M 2 : solid line - from NLC SR (||) with 
so = 2.4 GeV 2 , short-dashed line - from NLC SR without p'-meson contribution (so = 1.5 GeV 2 ), 
dashed line - from local SR of B&B§ (s = 1.5 GeV 2 ). 




Fig. 2: Extracted squared mass of the //-meson (in GeV 2 ): dashed line - sq = 2.9 GeV 2 , solid line - 
s = 2.6 GeV 2 , dotted line - s = 2.3 GeV 2 . 




Fig. 3: Longitudinal wave function of the /5-meson: solid line - from NL QCD SR., dotted line - 
asymptotic WF, dashed line - C&Z WF. 




Fig. 4: Longitudinal wave function of the //-meson: solid line - from NL QCD SR, dotted line - 
asymptotic WF. 



